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Abstract. Competition is a main tenet of economics, and the reason is that a 
perfectly competitive equilibrium is Pareto-efficient in the absence of externalities 
and public goods. Whether a product is selected in a market crucially relates to 
its competitiveness, but the selection in turn affects the landscape of competition. 
Such a feedback mechanism has been illustrated in a duopoly model by Lambert 
et ah, in which a buyer’s satisfaction is updated depending on the freshness of a 
purchased product. The probability for buyer n to select seller i is assumed to be 
Pn,i oc where Sn,i is the buyer’s satisfaction and T is an effective temperature to 

introduce stochasticity. If T decreases below a critical point Tc, the system undergoes a 
transition from a symmetric phase to an asymmetric one, in which only one of the two 
sellers is selected. In this work, we extend the model by incorporating a simple price 
system. By considering a greed factor g to control how the satisfaction depends on 
the price, we argue the existence of an oscillatory phase in addition to the symmetric 
and asymmetric ones in the (T, g) plane, and estimate the phase boundaries through 
mean-field approximations. The analytic results show that the market preserves the 
inherent symmetry between the sellers for lower T in the presence of the price system, 
which is confirmed by our numerical simulations. 
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1. Introduction 

Socio-economic systems have attracted attention of physicists due to their inherent 
dynamic complexity [iliaii. Simple interactions between social agents are known to 
give rise to non-trivial collective patterns like, e.g., segregation [ll|5l[6], opinion [ZllHlIH] 
or language dynamics [lO], and crowd behavior mm- Focusing more specihcally on 
economic aspects, the interaction of sellers and buyers in a market can be viewed as a 
dynamical process. For instance, taking space into account and including transportation 
costs in the agents’ utility leads to the well-known Hotelling model [131 E]) in which 
stores try to hnd the optimal location to maximize their proht. 

Here we consider a different problem, neglecting spatial aspects but taking into 
account a hnite lifetime of products, leading to a potentially rich dynamics. Buyers 
and sellers change their respective states upon every purchase: Namely, each buyer 
evaluates the sellers based on the purchased products, and each buyer also updates the 
list of products in stock. Under certain condition, their interaction may form a positive 
feedback loop in such a way that a seller, if selected by a buyer, becomes more likely to 
be selected in the future. Such a mechanism has been proposed and analyzed in detail 
by Lambert et al. [15] : They have considered two sellers dealing with perishable goods, 
so that a seller can replace existing products with fresh ones when buyers continuously 
make purchases from the seller. As a result of this positive feedback, the symmetry 
among the sellers gets broken spontaneously, leading to a virtual monopoly. Lambert 
et al. have analytically identihed a critical threshold for this phenomenon, when the 
control parameter is the degree of randomness in buyers’ choices. 

However, their duopoly model in Ref. [15] does not contain any price system, and it 
is an interesting question whether the market can restore the symmetry when the sellers 
are allowed to charge different prices. If the price declines as time goes by, and if buyers 
are sensitive to the price as well as freshness, it is indeed plausible that the market can 
better resist the tendency toward a monopoly. We thus extend the model of Ref. [15] 
by including a price mechanism, which dehnes the price of a product as a function of 
its freshness, and check to which extent this prescription stabilizes the market. 

In this paper, we show that the symmetry is recovered at a lower degree of 
randomness in buyers by including a price system coupled to the freshness of perishable 
goods. The threshold is estimated on a mean-held level in the sense that we neglect 
huctuations among buyers. We also hnd the possibility of a third phase, in which 
buyers seesaw between the two sellers. The period of this seesaw motion is estimated 
approximately by assuming homogeneity of the products in stock. These analytic results 
are consistent with a phase diagram obtained numerically. 

This work is organized as follow: We explain our model in Sec. |2j In Sec. 13.11 
we analyze its time evolution under the assumption that buyers remember the past for 
a long time. The threshold for symmetry breaking is given by examining stationary 
states in Sec. 13.21 Section l3R] explains the origin of the oscillatory phase and gives an 
approximate estimate of the period. We briehy discuss implications of our hndings in 
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2. Model 


Suppose that two sellers are competing to sell products of the same kind. Each seller 
has Np products in stock, and each of them has its own age r. As soon as a product 
is sold, it is replaced by a new product with r = 0 to conserve the number of products 
on the market all the time. As in Ref. [15], our assumption is that the products are 
perishable, so that the freshness changes with r in the following functional way: 

h(T) = e-^>^, ( 1 ) 


where Ti is the characteristic time scale for aging. It is a common marketing strategy 
to lower the prices of shopsoiled products. We may therefore assume that each seller 
can choose a price policy, parametrized by hp. This parameter means the characteristic 
freshness for a markdown. If he is low, the price will not drop down until the product 
becomes very old. Specihcally, we set 

x{h{T)) = x{t) = 1 — exp 


( 2 ) 


which is a decreasing function of r, bounded between zero and 1 — e~^^ . On the other 
hand, we have Na buyers. Let n G {!,..., A^^q} be an index to denote each buyer. If 
buyer n’s satisfaction from the Ah seller is denoted by the probability to choose 
this seller is given as 


Pn,i 


V ^pSn,il'^ 


( 3 ) 


where T is an effective temperature to control the degree of stochasticity in making the 
decision, and is a normalization factor. This formalism is often 

called the logistic selection model or logit rule [16], and it bears some analogy with the 
Glauber rate in physical systems, if we assume that energy plays the role of negative 
satisfaction. Suppose that the buyer n chooses the Ah seller to purchase a product 
of freshness hn,i for price Xn^i at time t. This product updates his or her satisfaction 
from this seller and affects the buyer’s next purchase at t' > t. We assume that the 
satisfaction at t' can be described as 

SnA^') = o^SnAt) + (!-«) {g (1 - xA) + (1 - 9)Ki,i} ) ( 4 ) 

where a characterizes the buyer’s memory and g is a. greed factor to determine the 
sensitivity to the price. The superscript t means that the price and freshness are 
measured at time t. 


3. Analysis 

3.1. Time evolution of probabilities 

Let us assume that the buyers can be represented by a single average buyer. In other 
words, we will neglect fluctuations in the buyers’ satisfaction as well as their behavior 
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to drop the buyer index n. Then, the probability to choose the ith seller becomes 
Pi = with Z = ^At time t, our average buyer pays x* to purchase a 

product of freshness hj from the ith seller. According to Eq. (j4]), the buyer’s satisfaction 
evolves as 

Si{t + To) = aSi{t) + (!-«) {g (l - x-) + (1 - g)hl} , (5) 

where tq denotes the average time interval to make a purchase. If the memory decays 
very slowly, i.e., 0<1 — a<Cl, we may regard S', as a smooth function of t to obtain 

• 1 — Ot — 

Si = - [g (1 - Xi) + (1 - g)hi - S'*] , (6) 

'To 

where the dot means the time derivative. Note also that we have suppressed the 
superscript t, because Eq. ([6]) is not a difference equation between t and t + tq, but 
a differential equation at time t. The time evolution of Pi can also be expressed as 



= 4p. (Si - {S}) , (7) 

where {X) = J2jPj^j means the average over the sellers. Equation ([7]) can be 
understood as the replicator equation [niiisiiis], where Eq. ([2]) plays the role of htness. 
By substituting Eq. (j6]) into Eq. ([7]), we hnd 

1 — _ _ 

Pi = — 7 ^Pi { [^(1 - Xi) + (1 - g)hi - Si] - {g{l - x) + (1 - g)h - S)} 

1 — (y. — 

= — 7 ^Pi {[gy - Xi) + {I - g)hi - Tin Pi] 

XqJ- 

-{g{l-x) + {l-g)h-Thip)] , (8) 

where we have used Si = T {In pi + InZ). 

3.2. Stationary state 

We hrst consider a stationary state in which all pfs are constant. The average fraction 
of products purchased from the ith seller during a time interval (t, t + dt) is expressed 
as 

Xi{t)dt = ^^f^dt, (9) 

NpTo 

because the average number of buyers visiting i per unit time is NaPi/xQ. The 
distribution of product ages r at the ith seller is denoted as ffr, t) with a normalization 
condition dT(j)i{T,t) = 1. The buyer randomly chooses a product, regardless of its 
age. The time evolution of cf)i{T,t) is thus written as 

(f)i{T + dt, t + dt) = 0i(r, t) — [\i{t)dt] (j)i{T, t). 


( 10 ) 
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By expanding Eq. fIlOp to the linear order of dt, we find that 

-Ai(t)0i(r,t). (11) 

This equation has a stationary solution ffir) = Xie~^*'^, from which hi for the 
stationary state is obtained as hi = dT</)f(T)h(T) = 1/ [1 + (AjTi)”^]. If we introduce 




R = 


— _ 


NpTQ 


= XiTi/pi, this result can also be written as 


hi — { 1 + 


Rpi 


-1 


( 12 ) 


In addition, the value of Xi = dT(l)f{T)x{T) can be computed from the stationary 
distribution as follows: 

1 — exp 


Xi = 


drXiC 


— XiT 


hr 


(13) 


Let us introduce another variable y = e The change of variables then leads to 

Xi = l- PiRh^f^'j [piR, hf^) , (14) 

where 'y{s,x) is the lower incomplete gamma function 'y{s,x) = dye~^y^~^. 

A possible stationary state is such that each seller is chosen with equal probability, 
i.e.. Pi = pq = 1/2. We will check the stability of such a state under small perturbation 
in the following way: First, let us rearrange Eq. ([8]) as 


X Of _ _ 

Pi = [g [Pi(l - Xi) -Pi{l- x)] + (1 - ^) [pihi - Pi (h)] 

- T[pi\npi-pi{\np)]} . (15) 

By adding small perturbation, we change the probabilities to Pi = po + Ci with = 0- 

The expressions inside the square brackets on the right-hand side (RHS) of Eq. flTHl) can 
then be expanded as 

Pi{l-Xi) -pi{l-x) = R{R,po,hc)ei +0{e^i), 

hiPo _ , 2 ^ 


(l + Rpo)^ 

= ei + 0(ef), 


0(ef), 


with 


Pihi - Pi (h) 

Pi In Pi -pi (Inp) 

r { R , Po, Ac) = T {pi)(p„ + e)iiA<>’»+‘>“7 [(po + €)R, hp ] } 


( 16 ) 

( 17 ) 

( 18 ) 


(19) 


e=0 


We note that the zeroth order of e* does not exist, so that Eq. flTHl) can be rewritten to 
the linear order of e,- as 


a 


roT 


gR{R,Po,hc) + (1-5')- 


Rpo 


T ^ ei. 


( 20 ) 


(1 + Rpo)^ 

The unperturbed stationary state is stable when the prefactor of e* on the RHS is 
negative. In other words, we can observe a symmetric phase, in which pi = 1/2 for 
every i, provided that 

Rpo 


T > Tc{g) = gR{R,po, K) + {1 - g) 


(l + Rpo)^' 


( 21 ) 
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for given R and q. One should note that multiple stable states can coexist in some range 
of R [15]. If R is small enough, however, we may say that the symmetric state is the 
only possibility when T > Tc. The question is what happens when T becomes lower 
than Tc- When the greed factor g approaches zero, the above analysis reduces to the 
result of Ref. na, in which an asymmetric phase emerges below Tc{g = 0). The origin 
of the asymmetric phase can be explained as follows: If the buyer’s choice gets slightly 
biased against a certain seller by chance, the products of this seller get older, which 
in turn lowers the probability to choose this seller, forming a positive feedback loop. 
In the end, only a single seller occupies the whole market, which is a stable stationary 
state at low T. Even if g has a small hnite value, it is reasonable to expect the same 
spontaneous symmetry-breaking phenomenon. When g gets high enough, however, the 
price system comes into play in a nontrivial way as will be explored below. 

3.3. Oscillatory phase 
Let us rewrite Eq. (|6|) as 

S^ = ^^[-S, + Qi{T)], ( 22 ) 

'To 

by dehning Qj(r) = g[l — Xi{T)] + {l — g)hi{T). It implies that we may expect Si ~ Qiir) 
when the system is evolving slowly. Suppose that only one seller, say, i = 1, is 
selected due to the symmetry-breaking mechanism explained above. This seller is always 
equipped with new products, so the satisfaction from this seller will be Si ~ Qi(0). On 
the other hand, r increases for the other seller as time t goes by, so 5*2 (^) will follow 
the trajectory of Q 2 (t). Figure [Dshows two possible cases of Qi(r), i.e., with small and 
large greed factors, respectively. When g is small, Q 2 {t > 0) stays below Qi(0), so the 
asymmetric phase remains stable [Fig. [11(a)]. On the other hand, if g is large enough, 
the low price compensates for the freshness reduction, and Qiir) can actually be greater 
than (5j(0) when r exceeds a certain threshold [Fig. [T](b)]. Now, the second seller begins 
to be selected, and the situation is less advantageous to the hrst seller: The value of Qi 
decreases, because some products are not sold, which lowers the freshness further. On 
the other hand, Q 2 does not decrease, because old products make positive contributions 
to Qi due to their low prices while refreshed products also attract buyers. The result is 
that the second seller is preferred until the prices of the other seller become low enough. 
To sum up, they take turns in playing the role of the champion and the challenger, 
which is the origin of the oscillating behavior. 

The above argument leads to some quantitative predictions: Assume that the 
behavior can be approximated as that of the asymmetric phase on a short time scale. We 
therefore suppose that p* = 1 so that only seller i is being selected. In this case, h and x 
are derived in Eqs. (IT^ and (fldjl . respectively, whereby we obtain an explicit expression 
for Qi at Pi = 1. This is denoted as Qq and represented as dotted horizontal lines in 
Fig. m Note that Qq ~ Qi(r = 0), because, by assumption, this sellers’ products are 
being replaced with fresh ones. However, it takes some hnite time to sell the products. 
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Figure 1. Q{t) for (a) g = 0.3 and (b) g = 0.6. The other parameter values are 
as given in Sec. IH The dotted horizontal lines represent Qo, and the existence of the 
second crossing point between Q(t) and Qq distinguishes the oscillatory phase from 
the asymmetric phase (see text). 

SO Qq is slightly lower than Qiir = 0). As a result, Qiir) has a trivial crossing at 
r ~ 0{NpTo/Na) (Fig. [T]). For old products to be as competitive as fresh ones, Qj(r) 
should have another crossing with Qq at some t = t* ^ 0{NpTo/Na). The shape of the 
curves in Fig. [Usuggests that this condition can be rephrased as Qj(r -P- oo) = g > Qq. 
It is clear that Qi(r) approaches g as r grows, so that the critical value gc is a solution 
of the following equation: 

^ [l - Rhf-f {R, 1/ he)] + = 9, (23) 

where the left-hand side expresses Qq. We may expect oscillatory behavior when g 
exceeds gc- At the same time, we stress that this calculation assumes deterministic 
behavior when it comes to buyers, so that the prediction will be the most accurate at 
T = 0. Furthermore, the crossing point r* between Qj(r) and Qq provides an estimate 
of the half-period of oscillation. It is obtained by numerically solving the following 
equation: 

^7exp + (1 - = Qq. (24) 

4. Numerical results 

Our numerical simulation shows the existence of the symmetric, oscillatory, and 
asymmetric phases as predicted above (Fig. [2]). In this simulation, we set the number 
of products at each seller as Ap = 5 x 10^ and the number of buyers as Na = 10^. The 
typical time interval of a purchase is given as tq = 10“^, and freshness is assumed to 
decay with a characteristic time scale ri = 20. The price policy is parametrized by a 
freshness threshold he = 5 x 10“^, and the memory factor is set as a = 0.99. With ah 
these parameters hxed, we change g and T to locate the phase boundaries. Recall that 
the boundary of the symmetric phase is predicted by Eq. fl?T]) . whereas the boundary 
between the asymmetric and oscillatory phases is predicted to be = ge, obtained by 
solving Eq. (|23|) . 

Let us dehne some numerical quantities to distinguish the three phases. First, the 
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Figure 2. Time evolution of the probability to choose each seller (red or black line), 
started from an initial condition pi = 1/2. As the greed factor g varies at a sufficiently 
low temperature T = 2 x 10“^, we observe (a) asymmetric, (b) oscillatory, and (c) 
symmetric phases, respectively, (d) With a relatively high temperature (T = 6 x 10“^), 
we can also find asymmetric oscillation. 



Figure 3. Phase diagram from numerical simulation (see text). For comparison, 
the phase boundaries predicted in Eqs. m and (1^ are represented by the solid 
lines. The labels S, A, and O mean the symmetric (m 0), asymmetric (m > 0), 
and oscillatory (m < 0) phases, respectively. In the region labeled as A’, we observe 
asymmetric oscillation [Fig. [2Kd)]. 
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g 

Figure 4. r* in Eq. (l24l) and numerically estimated half-periods, for different values 
of g at T = 10“^. The other parameters are the same as used in Sec. S] Equation (l24l) 
gives a good prediction for the half-periods when g is large enough, but the estimate 
becomes inaccurate as g approaches gc ~ 0.44, which is represented by the dotted 
vertical line. 

symmetry breaking can be detected by the following: 

mA=\{pi{t)-P2{t))t\, (25) 

where pi{t) = N~^ J2nPn,i is the probability to choose seller i at time t and means 
the average over time. Clearly, this quantity will be nonzero only in the asymmetric 
phase. To make a distinction between the symmetric and oscillatory phases, we also 
measure the following: 

mo = (bi(t) -P2it)\)t - \{piit) -P2{t))t\, (26) 

which will be nonzero only in the oscillatory phase. If we combine these two into a 
single parameter m = — mo, its value will be positive (negative) in the asymmetric 

(oscillatory) phase, and close to zero in the symmetric phase. Figure [3] shows the 
results: The solid lines are the phase boundaries estimated by Eqs. and fl2^ . The 
boundary between the symmetric and asymmetric phases shows an excellent agreement, 
especially for g < 0.8. The division between the oscillatory and asymmetric phases is 
also consistent with the estimation of gc when T is low. The region A’ is characterized 
by asymmetric oscillation, which has not been captured by our analysis. Such a phase 
can exist at sufficiently high T, because the unfavored seller sells its old products too 
early due to buyers’ stochastic choices. 

Finally, we compare r* with the half-period in the oscillatory phase in Fig. IH We 
choose very low temperature, T = 10“^, because the estimation of r* can easily be 
disturbed by noise. As shown in Fig. 01 the agreement is striking for large g. The 
failure at small g is related to one of our fundamental assumptions: Recall that we have 
assumed that all the buyers are homogeneous by using a single average buyer in Sec. 13.11 
When g is close to gc ~ 0.44, this assumption breaks down and the buyers have different 
probabilities to choose the sellers. This inhomogeneity dramatically reduces the period 
of oscillation. 







Symmetry restoration by pricing in a duopoly of perishable goods 10 

5. Discussion 

The phase diagram in Fig. |3]shows that the price system indeed restores the symmetry in 
a broader range of T as compared to the case withont price formally obtained as 5 ^ = 0. 
However, the asymmetric phase can exhibit oscillations when bnyers are sensitive to the 
price. One might point ont that onr model bears similarity to the Bertrand dnopoly, 
in which two sellers sell homogeneons prodncts withont any possibility of collusion [ 20 ] . 
The sellers in the Bertrand duopoly compete by setting prices simultaneously, and the 
one with a lower price occupies the whole market. If the two sellers charge the same 
price, they will evenly share the market. The pure-strategy Nash equilibrium turns out 
to be such that both of the sellers set the price equal to marginal cost of production, 
which means that they earn nothing from their business [21]. In the Bertrand duopoly, 
therefore, the existence of two sellers is already enough to actualize the lowest possible 
price. However, one should note its underlying assumption that a single seller can cover 
the whole market. If this is not the case, e.g., due to limited capacity of production, 
we encounter the Edgeworth paradox, which means that this game has no pure-strategy 
Nash equilibrium, hence no stationary price. For example, one of the two hrms can make 
prohts by deviating from the marginal cost, because the other hrm alone cannot meet 
the demand. Although a mixed strategy can constitute a Nash equilibrium [2^ |23], it 
is hardly feasible in practice and one would instead observe the Edgeworth price cycle 
(Reference [20] also points out that the capacity constraint may not be a necessary 
condition for the cycle). The Edgeworth paradox suggests that the price mechanism 
can introduce instabilities. However, the cycle is actually a pseudo-dynamic process in 
the sense that it originates from an equilibrium concept. In contrast, our oscillatory 
phase is born out of dynamical rules without any consideration of strategic equilibrium. 

To consider some strategic aspects in our model, let us suppose that the variables 
T and g characterizing buyers are just given parameters when viewed from the sellers. 
The sellers can instead decide he to control the price policy, and the symmetry between 
the sellers suggests that they will end up with the same he- The critical temperature 
Tc in Eq. (|2T]) depends on he unless g = t). For example, if he —)■ 0, the price will 
become insensitive to freshness. It implies that iF{R,po, he) vanishes in Eq. fl2l|) so that 
the phase boundary of the symmetric phase in Fig. |3]will converge to a line connecting 
(Tc(0),0) and (0,1) on the (T,g) plane. In a price war, the sellers will increase he to 
lower the average price, but the possible range of he must be bounded due to the cost 
of production. We may regard our he in the previous sections as the highest possible 
one determined by such a competitive process. Differently from the Bertrand duopoly, 
however, the sellers do not always divide the market half and half, when the mechanism 
of Lambert et ah is at work [15]: The market of perishable goods tends to become 
monopolistic, especially when the purchasing behavior is deterministic with low T. Our 
point is that this tendency is nevertheless weakened by the price sensitivity, either by 
extending the symmetric phase, or by restoring the symmetry over a long period in the 
oscillatory phase. 
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6. Summary 

To summarize, we have incorporated a simple price system coupled with freshness into 
the duopoly model suggested by Lambert et ah In addition to the effective temperature 
T to control the randomness in purchasing behavior, we have introduced the greed 
parameter g to determine the sensitivity of satisfaction to the price. We have identihed 
the symmetric, asymmetric, and oscillatory phases and estimated their boundaries in 
the (T, g) plane. Our numerical simulations show nice agreements with our analytic 
results. Based on our analysis, we conclude that the price system resists the tendency 
to a monopoly: On one hand, it lowers the critical temperature Tc below which the 
symmetric phase becomes unstable. On the other hand, the market with high g can 
keep oscillating without settling on a single seller, preserving the symmetry in a time- 
averaged sense. 
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